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QCD SUM RULES STUDY OF THE J^'^ = l~ CHARMONIUM Y MESONS 
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C.P. 66318, 05389-970 Sao Paulo, SP, Brazil 

We use QCD sum rules to test the nature of the recently observed mesons 1^(4260), y(4350) 

and y(4660), assumed to be exotic four-quark (ccqq) or (cess) states with J = 1 . We work 

at leading order in as , consider the contributions of higher dimension condensates and keep terms 

which are linear in the strange quark mass tUs. We find for the (cess) state a mass my ~ (4.65±0.10) 

00 ' GeV which is compatible with the experimental candidate y(4660), while for the (ccqq) state we 

c7^ I find a mass my = (4.49 ±0.11) GeV, which is bigger than the mass of the experimental candidate 

^^ . y(4350). With the tetraquark structure we are working we can not explain the F(4260) as a 

C^ ' tetraquark state. We also consider molecular DsqDI and DqD* states. For the DsoD* molecular 

5^ ' state we get m£)_,g£|. — (4.42 ±0.10) GeV which is consistent, considering the errors, with the mass 

&*-(■ of the meson y(4350) and for the DqD* molecular state we get mjj^f), = (4.27 ± 0.10) GeV in 

■^T ' excelent agreement with the mass of the meson Y(4260). 

o 
m 

I ■ I. INTRODUCTION 

^: 

(~| I Recent experimental observations of new particles at i3-factories, have revitalized the field of hadron 

spectroscopy. While some of these observed states, like the X(3940) [ll], can be understood as cc states 
[31, these observations suggest that the spectrum of the charmonium states is much more rich than 
suggested by the quark-antiquark model and may include tetraquark or molecular states. In particular, 
the X(3872) Q , with quantum numbers J^*^ = 1++, was the first observed state that did not fit easily 
the charmonium spectrum. It also presents a strong isospin violating decay that disfavors a cc assignment. 
00 ■ The Z+(4430), recently observed in the Z'^ -^ ip'ir^ decay mode [4], is the most intriguing one since, 

■^ I being a charged state, it can not be a pure cc state. 

t;;}-' '. Besides the X(^872) and Z+(4430) mesons, at least three peaks of the J^'^ = 1 — family: the r(4260) 

O ; [I, the y(4325) [| or/and r(4360)J7i (that here we call r(4350)) and the y(4660) 0] can not be easily 

00 , fited in the charmonium spectrum Q- Since the only observed decay channels for these states are those 

^^ ' containing the J/tp (for F (4260)) or ■0' (for the others) plus a pair of pions, in ref. [2] this was considered 

as an indication that these mesons contain a particular charmonium resonance, J/tjj or ip', that stays 
intact inside a more complex hadronic structure. 

A critical information for understanding the structure of these states is wether the pion pair comes from 
^ ■ a resonance state. From the di-pion invariant mass spectra shown in ref. [^ there is some indication that 

only the y(4660) has a well defined intermediate state consistent with /o(980) ¥|. Due to this fact and 
the proximity of the mass of the ?/'' — /o(980) system with the mass of the y(4660) state, in ref. 1Q\, the 
y (4660) was considered as a /o(980) V'' bound state. The y(4660) was also suggested to be a baryonium 
state [11| and a canonical 5 '^Si cc state J12l ). 

In the case of y(4260), in ref. [l3| it was considered as a sc-scalar-diquark sc-scalar-antidiquark in 
a P-wave state. In a naive estimate, the mass of a sc-scalar-diquark would be approximately equal to 
the mass of the D^ meson. Since a unit of angular momentum can be estimated as the mass difference 



between the nucleon and the lowest-lying odd parity excited state, and it is about 600 MeV [lj|, one 
would expect the mass of a sc-scalar-diquark sc-scalar-antidiquark in a P-wave state about 4540 McV, 
almost 300 MeV above the y(4260) mass. 

Here we use the QCD sum rules (QCDSR) [1^, [la, Il3| to study the two-point function for a diquark- 
antidiquark tetraquark state with a symmetric spin distribution: [cs] 5=1 [cs]s=o + [cs]s=o[cs]s=i, to 
see if any of these new Y mesons can be described by such a current. In previous calculations, the 
QCDSR approach was used to study the AT (3872) considered as a diquark-antidiquark state [l^], and the 
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Z+(4430) meson, considered as a D*Di molecular state [Ta|- In both cases a very good agreement with 
the experimental mass was obtained. 

II. THE TWO-POINT CORRELATOR 

The lowest-dimension interpolating operator for describing a J^^ = 1 state with the symmetric 
spin distribution: [cs]s=o[cs]s=i + [cs] 5=1 [cs] 5=0 is given by: 

Jm = -^^?|^[(SaC'75C6)(srf7^75Ccf ) + (srC757^C6)(sd75Ccf )] , (1) 

where a, b, c, ... are color indices and C is the charge conjugation matrix. 
The two-point correlation function is given by: 

n^,{q) ^t fd^x e''--(0|r[j4x)jt(0)]|0) = -Il{q'){g^,q' - q,q,), (2) 



from where we get 

U^^iq) = -3q'll{q'), (3) 

The calculation of the phenomenological side at the hadron level proceeds by writing a dispersion 
relation for the invariant function in Eq. Q: 

nP''™(g2) = - fds , ^^''^ . + ■■■ , (4) 

J q^ - s + ie 

where p is the spectral density and the dots represent subtraction terms. The spectral density is described, 
as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum representing 
higher mass states: 

p{s) = A2<5(s-m2.)+p-"*(5), (5) 

where A is proportional to the meson decay constant fy, which parametrizes the coupling of the vector 
meson, Y , to the current j^: 

(OlJ^li") = /yTO^e^ = Xmyef, . (6) 

For simplicity, it is assumed that the continuum contribution to the spectral density, /9'^°"*(s) in Eq. ([5]), 
vanishes bellow a certain continuum threshold Sq. Above this threshold, it is assumed to be given by the 
result obtained with the OPE. Therefore, one uses the ansatz [20| 

p— *(s) = pO^^(s)e(s-so), (7) 

with 

pOPE^^-) ^ l/m[nOP^(s)] . (8) 

TT 

In the OPE side, we work at leading order in a^ and consider the contributions of condensates up to 
dimension six. We keep the term which is linear in the strange-quark mass mg. We use the momentum 
space expression for the charm quark propagator, while the light-quark part of the correlation function is 
calculated in the coordinate-space. The correlation function, n(g^), in the OPE side can also be written 
in terms of a dispersion relation: 

U^P^iq')^ ds^- ^. (9) 

J Ami ^ ~ 9 

Therefore, after making an inverse-Laplace (or Borel) transform on both sides, and transferring the 
continuum contribution to the OPE side, the sum rule for the vector meson Y can be written as 

JAml 



where 
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where the integration limits are given by amin = (1 ^ \/^ — 4to2/s)/2, amax = (1 + \/l — Am'^/s)/2 
and (3min — aml/{sa — ml). The contribution of dimension-six condensates {g^G^) is neglected, since 
it is suppressed by the loop factor 1/I67r2. In Eq. (|13p we have included, for completeness a part of the 
dimension-8 condensate, related with the mixed condensate-quark condensate contribution. We should 
note that a complete evaluation of these contributions require more involved analysis including a non- 
trivial choice of the factorization assumption basis [2l|. In Eq. P^ the n'""^'*^^(M2) term is treated 
separately because its imaginary part is proportional to delta functions and, therefore, can be easily 
integrated J22 |. 

It is very interesting to notice that the current in Eq. ([1]) does not get contribution from the quark 
condensates when nis — 0. This is very different from the OPE behavior obtained to the scalar-diquark 
axial-antidiquark current used for the X(3872) meson in ref. [l^, but very similar to the OPE behavior 
obtained for the axial double-charmed meson Tec, also described by a scalar-diquark axial-antidiquark 
current (2311. 



III. SUM RULE PREDICTIONS FOR my 



The values used for the quark masses and condensates are |17ll24|: m,c{m,c) = (1.23±0.05) GeV, m-s = 
(0.13±0.03) GeV, {qq) = -(0.23±0.03)3 GeV^ {ss) = 0.8{qq), {qgcr.Gq) = ml{qq) with ml = 0.8 GeV^ 
(52^2^ ^ 0.88 GeV^ 

We evaluate the sum rules in the Borel range 2.8 < M'^ < 4.6 GeV^, and in the sq range 5.0 < ^/sq < 5.2 
GeV. To fix the continuum threshold range we extract the mass from the sum rule, for a given sq, and 
accept such value of Sq if the obtained mass is around 0.5 GeV smaller than ^/sq. 

From Fig. [1] we see that we obtain a quite good OPE convergence for AP > 3.2 GeV2. Therefore, we 
fix the lower value of IVP in the sum rule window as Af^j„ — 3.2 GeV2. This figure also shows that 
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FIG. 1: The OPE convergence in the region 2.8 < M^ < 4.5 GeV for ^/s^ = 5.1 GeV. Perturbative contribution 
(soUd line with triangles), (ss) contribution (dashed-line) , {g^G^) contribution (dotted line), {sga.Gs) contribution 
(solid line with squares), (ss)^ contribution (dot-dashed hne), {sga.Gs){ss) contribution (sohd line with spheres) 
and total contribution (solid line). 
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FIG. 2: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for ^/so = 5.1 GeV. 



dimension-eight condensate contributions is very small. Since this is not a complete evaluation of the 
dimension-eight condensates contribution, it will be neglected in this calculation. 

In Fig. [2] we show the comparison between pole and continuum contributions for y/sQ = 5.1 GeV, 
and we see that for M^ < 4.05 GeV^, the pole contribution is bigger than the continuum contribution. 
Therefore, we fix M^ = 4.05 GeV^ as the upper limit of the Borel window for ^/sq = 5.1 GeV. The same 



analysis for the other values of the continuum threshold gives M < 3.8 GeV for ^/sq = 5.0 GeV and 
M^ < 4.2 GeV^ for y/so — 5.2 GeV. We then consider, for each value of sq, the range of A'P values from 
3.2 GeV^ until the one allowed by the pole dominance criteria given above. 

To extract the mass my we take the derivative of Eq. ^TU\i with respect to 1/M ^, and divide the result 
by Eq. (fTU|). In Fig. [31 we show the Y meson mass, for different values of ^/sq, in the relevant sum rule 
window, with the upper and lower validity limits indicated. From this figure we see that we get a very 
good Borel stability for my- 
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FIG. 3: The Y meson mass as a function of the sum rule parameter (M^) for different values of ^/sq. The arrows 
indicate the region allowed for the sum rules: the lower limit (cut below 3.2 GeV^) is given by OPE convergence 
requirement and the upper limit by the dominance of the QCD pole contribution. 



To check the dependence of our results with the value of the charm quark mass, we fix ^/s^ =5.1 GeV 
and vary the quark masses and the quark condensate in the range given above. We see that the results 
are not very sensitive to the changes in the values of the quark masses and condensate, the main source 
of uncertainty being the continuum threshold. Adding the errors in quadrature we finally get 



my = (4.65 ±0.10) GeV, 



(14) 



in excelent agreement with the mass of the y(4660) meson. Therefore we conclude that the meson 
y(4660) can be described by a tetraquark state with a spin configuration given by scalar and vector 
diquarks. 



IV. SUM RULE STUDY FOR F(4350) 

As pointed out in the introduction, from the di-pion invariant mass spectra shown in ref. [9|, there is 
some indication that only the y(4660) has a well defined intermediate state consistent with /o(980). The 
di-pion invariant mass spectra for the pions in the decay y(4350) — > ■(/''tt+tt" could be consistent with a 
broad intermediate state with a mass around 600 MeV. This could be a ct scalar meson [2F 
case, one would not expect strange quarks in the current in Eq. ([1]). 

Replacing the strange quarks in Eq. ([T]) by the generic light quark q and using ri 



and, in this 



7.0 MeV, we show, 
in FigUl the OPE convergence using y/so = 4.9 GeV 

Since now the quark condensate contribution is multiplyed by the light quark mass, its cotribution is 
now neglegible. From this figure we see that we get a good OPE convergence for A'P > 3.2 GeV^. 

Again the upper limits for Af ^ are obtained by imposing that the pole contribution should be bigger 
than the continuum contribution. These values are given in Table I for each value of y/sQ. 

Table I: Upper limits in the Borel window for the sum rule with nis — niq. 
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In Fig. [5] we show the relative continumm (solid line) versus pole (dashed line) contribution, using 
y/so = 4.9 GeV, from where we clearly see that the pole contribution is bigger than the continuum 
contribution for AP < 3.5 GeV^. 
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FIG. 4: The OPE convergence for the sum rule with rUs ~ ruq, using ^/sg — 4.9 GeV. The solid with trian- 
gles, dashed, dotted, solid with squares, dot-dashed and solid lines give, respectively, the perturbative, quark 
condensate, gluon condensate, mixed condensate, four-quark condensate and total contributions. 
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FIG. 5: Same as Fig. 2 for the sum rule with rus — niq and y^ = 4.9 GeV. 

In this case the stability for the ?7iy mass is not as good as in the previous case, but it is still acceptable, 
in the allowed sum rule window, as a function of M^ as can be seen by Fig. [6l 

Taking into account the variations on A/^, so, {qq) and rric in the regions indicated above we get: 

my = (4.49 ±0.11) GeV , (15) 

which is bigger than the y(4350) mass, but it is consistent with it considering the uncertainty. 



V. DsoDl MOLECULE 

As pointed out in ref. [26|, if the X(3872), and Z+(4430) are really molecular states, then many other 
molecules should exist. In particular a I?so(2317)Z)*(2110) molecule with J^*" — 1 , could also decay 
into '0'7r^7r~ with a dipion mass spectra consistent with /o(980). Therefore, in this section we consider 
a meson-meson current to check if this current could also describe the meson y(4660). A current with 
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FIG. 6: The light my meson mass as a function of the sum rule parameter (M^) for different values of ^/so- The 
arrows delimit the regions allowed for the sum rule. 



J — 1 and a symmetrical combination between the scalar and vector mesons is given by: 

jp = —^[{Sa^t,Ca){CbSt) + {Calt,Sa){SbCb)] ■ 



(16) 



With this current we get for the spectral density, results very similar to the ones in Eg. p^ . the only 
differences being the color factors. We get a factor 3/4 for the diagrams where there is no gluons being 
exchanged between the quark lines, and a factor 3/2 for the diagrams with gluons being exchanged 
between the quark lines: 
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With these changes, the mixed condensate is now, for the current in Eq. (jl6[) . more important than 
it was for the current Eq. ([T]), as can be seen in Fig. [T] For this current wc find that the appropriate 
continuum threshold range is 4.8 < ^/s^ < 5.0 GeV and that the OPE convergence is very good for 
M^ > 3.2 GeV . The upper limits for M^ are are given in Table II for each value of ^/sq. 
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FIG. 7: The OPE convergence for the sum rule for the DsoD* molecule with ^/so = 4.9 GeV. 



Table II: Upper limits in the Borel window for the sum rule for the DsoD* molecule. 
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FIG. 8: The DsoD^ molecule mass as a function of the sum rule parameter (M^) for different values of ^/sq. The 
arrows delimit the regions allowed for the sum rule. 



From Fig. [8] we see that the Borel stability, as a function of M^ is very good, in the allowed sum rule 
window. 

Taking into account the variations on M^, sq, (qq) nis and rric in the regions indicated above we get: 



iD^^a-i = (4.42 ±0.10) GeV 



(18) 



which is more in agreement with the F(4350) mass than with the F(4660) mass. It is important to 
mention that even if we use ^/sq = 5.1 GeV, which is the central value used for y(4660), we get 
''^DsoD- ~ 4.5 GeV, still in agreement, considering the errors, with the value in Eq. (fT8|l . Therefore, 



we have to conclude that the y(4660) meson is better explained with the diquark-antidiquark current in 
Eq. Il]) than with the molecular current in Eq. (flB]) . To conclude if we can associate this molecular state 
with the meson y(4350) we need a better understanding of the di-pion invariant mass spectra for the 
pions in the decay F(4350) —f iJ)'tt~^tt~ . From the spectra given in ref. p, it seems to us that the y(4350) 
is more consistent with a non-strange four-quark state than with a DsqD* molecular state. However, 
from the mass obtained from the sum rules we see that we can explain the 1"(4350) meson better as a 
DsqDI molecular state. 

It is also important to notice that the central mass obtained for the D^qD* molecule is very close to 
the I?so(2317)-D|(2110) threshold, indicating that this channel might be forbidden in the F(4350) decay. 



VI. Do-D* MOLECULE 

Similarly to what was done in section IV, we can consider a scalar-vector charmed mesons molecule: 
DqD* with J^*^ — 1~~. For this we have only to change the strange quarks in Eg. p^ by the generic 
light quark q. In this case we get a good OPE convergence for M^ > 3.2 GeV^ and the continuum 
threshold in the range 4.6 < ^/sq < 4.8 GeV. The upper limits for AP are are given in Table III for each 
value of ^/sq. 

Table III: Upper limits in the Borel window for the sum rule for the DqD* molecule. 
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FIG. 9: The DqD* molecule mass as a function of the sum rule parameter (M^) for different values of 
arrows delimit the regions allowed for the sum rule. 



The 



Again we get a very good Borel stability as a function of M^ , in the allowed sum rule window, as can 
be seen by Fig. [9l Considering the variations on Af^, sq, (qq) and nic in the regions indicated above we 
get: 

mjj^Q, = (4.27 ±0.10) GeV , (19) 

in excelent agreement with the mass of the meson y(4260). Again, to conclude if we can associate this 
molecular state with the meson y(4260) we need a better understanding of the di-pion invariant mass 
spectra for the pions in the decay V(4260) -^ J/i/^Tr+Tr". From the spectra given in ref. [9], it seems that 
the V(4260) is consistent with a non-strange molecular state DqD*. 

Regarding charmed non-strange scalar mesons, a broad enhancement (Dq) mass distribution has been 
observed [27j. Its mass is now compiled as f28| mug = 2352 ± 50 MeV. Therefore, the DqD* threshold 
is around 2360 MeV and is 100 MeV above the molecule mass in Eq. (fTT ^ 



10 
VII. CONCLUSIONS 

In conclucion, we have presented a QCDSR analysis of the two-point functions for tetraquark char- 
moniuni states with J^'-' — 1 , to study the charnionium Y mesons recently observed by BaBar and 
BELLE Collaborations. We have considered two kinds of currents: a diquark-antidiquark tetraquark state 
with a symmetric spin distribution [cg]s=i[c9]s=o + [cg]s=o[c9]s'=i (where the generic quark q can be ei- 
ther a light u or d quark, or a strange, s, quark), and a molecular state with a symmetrical combination 
between scalar and vector charmed mesons. 

Our findings indicate that the y(4660) meson can be very well described by a diquark-antidiquark 
(cs) (csVtetraquark state. This quark content is consistent with the di-pion invariant mass spectra shown 
in ref. Q, which shows that there is some indication that the y(4660) has a well defined di-pion inter- 
mediate state consistent with /o(980). 

In the case of the F(4350) meson, its mass can be reproduced, considering the errors, if we assume that 
it is a diquark-antidiquark {cq){cq) tetraquark state, or a DsqD* molecule. Since an indication from its 
quark content can be obtained from the di-pion invariant mass spectra in the decay Y{A350) —^ ip' n'^ Tr~ , 
we need a better understanding of it before we can reach a definite conclusion about its structure. 

We also found that the F(4260) meson can be very well described by a DqD* molecular state. 
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